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Abstract. A conceptually novel model - one in which stimulated (i.e. induced) resonance Rayleigh scattering causes ions of 
select species to become accelerated to very high terminal velocities in the solar winds of OB-type stars - is proposed to explain 
the general appearance of so-called P Cygni profiles that often dominate the vacuum ultraviolet (VUV) spectra of such stars. 
In the unit step of the proposed scattering process, a radially outgoing photon (frequency vi) that is arbitrarily blueshifted with 
respect to the rest frame frequency v„ of the resonance line associated with a P Cygni profile is absorbed from the illuminating 
star's continuum, a new photon (frequency v 2 ) that is very close to v„ and that propagates in the backwards direction (i.e. towards 
the star) is created, and the outwardly directed velocity v of the ion being accelerated is increased by an amount Av w 2/iv D /cm/, 
with all three events occurring simultaneously. A monochromatic wave at v 2 initially forms at some distance from the star, and 
becomes enormously amplified via the stimulated scattering process as it propagates radially inwards towards the star, all the 
while retaining a relatively high degree of monochromaticity. The high rate of stimulated scattering enables ions of the resonant 
species in the solar wind to become accelerated to terminal velocities as high as ~1000 km/sec. Since stimulated scattering 
processes are characterized by pump power thresholds, the model readily explains why only select species are accelerated to 
very high velocities in a given star's solar wind, and also why a dramatic P Cygni profile for a given ion species can often 
discontinuously be present or absent when spectra of stars varying only slightly in spectral type are compared. 
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1. Introduction 

So-called P Cygni profiles of select ion resonances (e.g. C IV 
/U1548, 1551; Si IV AA1394, 1403; NV ^1239, 1243 - see 
Fig. 1) are often the dominant features in the vacuum ultra- 
violet (VUV) spectra of O-type (Walborn et aZ. Il985t and B- 
type (Walborn et al. 1995) stars, especially giants and super- 
giants. These features were originally discovered in the late 
nineteen-sixties by Morton (Morton 1967) and his colleagues 
(Morton et al. 11969b with the use of rocket-based instrumen- 
tation, and have since been studied extensively by many as- 
tronomers utilizing several orbiting satellites that were success- 
fully launched in the years that followed. Especially produc- 
tive in obtaining P Cygni spectral data were the fully dedicated 
Copernicus and International Ultraviolet Explorer {IUE) satel- 
lites, and the Space Telescope Imaging Spectrometer (STIS) 
aboard the Hubble Space Telescope. 

In a typical P Cygni profile, the star's continuum level in 
a relatively broad spectral region that is blueshifted with re- 



spect to a given ion resonance line v is heavily absorbed. 
(Frequently the ion resonances are doublets, as is the case with 
the three prominent P Cygni profiles shown in Fig. [J) The 
measured frequency displacement from v a of the short wave- 
length limit vb of this absorbed region is normally used by as- 
tronomers to determine the maximum velocity v max of the cor- 
responding ions in the solar wind via the standard relationship 
vb - v„ = v„v m a X /c. From such Doppler-shift-based analyses 
of P Cygni profiles, astronomers have inferred that the corre- 
sponding ions are strongly accelerated radially away from the 
stars, reaching in some cases maximum velocities as high as 
two or three thousand km/sec. For a hot star, the maximum 
velocity deduced in this manner is frequently greater than the 
calculated escape velocity V e = (2GM/R) 1 ^ 2 at the radius R of 
the photosphere, implying an actual ejection of matter from the 
star into interstellar space. 

Included in Morton J1976I are special Copernicus VUV 
spectral scans of the 04 If supergiant £ Pup, recorded at very 
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Fig.l. VUV spectrum of the 09 supergiant HD 210809 
recorded with the Space Telescope Imaging Spectrometer 
(STIS) aboard the Hubble Space Telescope. (Spectrum down- 
loaded from the MAST Scrapbook interactive web site.) The 
three P Cygni profiles dominating the spectrum are (from left 
to right) due to N V, Si IV, and C IV ions, respectively. The ab- 
sorption at 1216 A is that of H-atom Lyman alpha. A multitude 
of narrow interstellar absorption lines are present. Wider stellar 
photospheric absorptions can also be seen. 



high spectral resolution (0.051 A FWHM), and with an excep- 
tionally high signal-to-noise ratio. For this very hot star, which 
has R/R Q w 20.3 and M/M Q * 100, it was calculated in Morton 
(1976) that V e =1370 + 160 km/sec. Between 920 A and 1750 
A, 13 different ion species were found to display P Cygni pro- 
files in £ Pup, many of them appearing very strong. For ex- 
ample, the ion C III was observed to have a strong P Cygni 
profile associated with its resonance line at sa 977 A. The blue 
edge of the heavily absorbed region appeared shifted from 977 
A by a: 2710 km/sec, prompting the conclusion that was made 
in Morton (1976) that C III ions in the solar wind of £ Pup 
are continually being ejected into space. However, as will be 
shortly pointed out, according to the nonlinear P Cygni model 
here being proposed, use of the relationship vb - v = v v max /c 
to determine the maximum ion velocity v max attained in a star's 
solar wind results in an overestimation of this quantity by al- 
most exactly a factor 2. 

As noted in Lamers & Cassinelli (1999), current stellar 
wind theories fall into three broad classes: radiative (i.e. line- 
driven) models, coronal models, and hybrid models. In radia- 
tive models, transfer of photon momentum to the gas is as- 
sumed to occur through the opacity of the many strong VUV 
resonance lines that are present. Increased acceleration is be- 
lieved to result from the progressive Doppler shifting of the line 
opacity into the unattenuated photospheric radiation field. The 
stellar wind in existing radiative models is thus assumed to be 
driven by a purely linear effect, commonly known as radiation 
pressure. 

In the model here being proposed, a nonlinear photomecha- 
nism drives the acceleration of the fast moving ions in the solar 
wind, simultaneously producing the P Cygni spectral profiles 
one sees in the line-of-sight to the star. Perhaps the most at- 
tractive feature of such a nonlinear mechanism would be that 
it would possess a definite pump power threshold. Only those 
ion species for which this threshold is reached would be accel- 



erated to very high velocities in a given star's solar wind. Such 
selectivity is very hard to explain with linear radiative models 
for solar winds. 

In the present paper the focus is placed upon the basic 
physics of the proposed nonlinear model. No attempt is made 
to account in detail for fine-grained features appearing in P 
Cygni profiles of specific stars. The paper is organized as fol- 
lows. In Sect. 2, C II /L11334, 1336 ion resonance spectra 
recorded in four different stars of roughly the same spectral 
type are compared to demonstrate the abrupt manner in which 
P Cygni profiles often appear in star spectra, thus strongly sug- 
gesting that a threshold effect of some kind must be involved. 
In Sect. 3, the proposed nonlinear scattering mechanism is out- 
lined. Undoubtedly of some significance is the fact that this 
mechanism is easily shown to be consistent with the conserva- 
tion of both energy and momentum. Simple consideration of 
the manner in which spontaneous resonance Rayleigh scatter- 
ing (i.e. elastic scattering) would occur in the rest frame of an 
ion being accelerated in a star's solar wind directly leads to 
identification of the proposed nonlinear scattering mechanism 
with stimulated (i.e. induced) resonance Rayleigh scattering. In 
the unit step of this stimulated process, three distinct events si- 
multaneously occur. An outwardly propagating photon from 
the star's continuum that is blueshifted with respect to the rest 
frame frequency v of a strong ion resonance line becomes ab- 
sorbed by the ion. (0 A photon at exactly v that propagates 
radially backwards (i.e. towards the star) is emitted by the ion, 
contributing to the intensity of a backward propagating light 
wave at v„ that is already present at the position of the ion. (For 
analytical purposes, this wave is termed the V2 wave, but its fre- 
quency is k v .) The outwardly directed radial velocity v of 
the ion is increased by Av « 2hv /cmj, irrespective of the value 
of v. 

Since, at any point in the solar wind, the transition proba- 
bility for the unit step to occur is proportional to the product 
of the outwardly propagating continuum flux and the inwardly 
propagating light beam intensity 1%, the likely occurrence of a 
stimulated scattering regime is therefore here expected, with 
the V2 wave continuously gaining in intensity as it propagates 
inwardly towards the star. During the entire time it undergoes 
amplification, the V2 wave remains relatively monochromatic, a 
result of the unit step in the scattering process becoming stim- 
ulated by the same wave. However, as the V2 wave propagates 
inwardly towards the star, the frequency of the pump light that 
effectively drives the stimulated scattering process continually 
changes. In this way, continuum light from the star spanning 
a broad spectral range can contribute to the intensity of the 
monochromatic V2 wave as the latter impinges upon the pho- 
tosphere of the star. 

In Sect. 4, equations that should describe the stimulated 
scattering process are briefly discussed. However, efforts to 
provide numerical solutions for these equations have not yet 
been attempted. 

In the next section of the present paper (Sect. 5), two types 
of spectroscopic evidence that could help to substantiate the 
nonlinear model are discussed - fluorescence spectra and two- 
photon absorption bands. Most P Cygni profiles are observed 
to possess a prominent fluorescence component. Via simple ex- 
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tension of the stimulated Rayleigh scattering model used to 
explain P Cygni absorption, one can also comprehend how P 
Cygni fluorescence is excited. By considering the effect of hav- 
ing the absorption and the fluorescence excitation mechanisms 
operate in tandem, one can then easily account for the observed 
fact that the strongest apparent absorption in a P Cygni profile 
occurs at frequencies considerably blueshifted from v„. 

In Sect. 5 it is also suggested that one might seek to iden- 
tify Doppler broadened absorption bands representing reso- 
nantly enhanced two-photon absorption of outwardly propagat- 
ing continuum light from a star by various ion species present 
in its photosphere, or just outside. Such two-photon absorption 
could in principle be induced by the powerful, inwardly propa- 
gating, monochromatic beams at the P Cygni ion rest frame fre- 
quencies v a , provided that near resonances exist between some 
of the frequencies v and some transitions of the two-photon- 
absorbing ions. Details of a search for such near resonances in 
the case of the P Cygni star f Pup are described. However, with 
one possible intriguing exception, this search did not result in 
the clear identification of any two-photon absorption bands in 
this star. 

In Sect. 6 nonlinear and linear P Cygni mechanisms are 
briefly compared. 



2. P Cygni profiles in stars of almost similar 
spectral types 

In Fig. |5J VUV spectra recorded with the IUE satellite of four 
B-type supergiants are displayed over a limited wavelength 
range that includes the strongly allowed C II ion resonance 
lines at 1334.5 A and 1335.7 A. Each star shown is of a dif- 
ferent spectral type, but the difference in temperatures between 
stars in adjacent spectra is relatively small. In the figure, star 
temperatures decrease in going from top to bottom. 

It is seen that in (a) and (b) the C II doublet appears nor- 
mally absorbing. Whether the absorbing ions in these two spec- 
tra are located in the stars' photospheres, or just outside, is 
somewhat difficult to judge. It may also be that the lines, at least 
in part, represent interstellar absorptions. Whatever the case, 
a significant spectral change is seen to occur in (c), namely, 
a P Cygni profile, with its characteristic blueshifted region of 
strong absorption, dramatically appears. From all four spectra, 
it is clear that this intense new absorption is not simply a broad 
feature of the stellar photosphere. To explain its sudden appear- 
ance via theories based upon the effects of linear radiation pres- 
sure also seems excessively difficult. Linear photoexcitation is 
probably entirely responsible for the comparatively narrow C II 
absorption lines seen in (a) and (b). These lines are still seen to 
be present in (c) and (d), and linear absorption may still be the 
agent that produces them. However, it is quite evident from the 
spectral sequence shown in Fig. |3 that a dramatic new mech- 
anism significantly increasing the degree to which continuum 
light emitted by a star is attenuated can sometimes abruptly 
become activated. In the present paper it is proposed that this 
mechanism is stimulated Rayleigh scattering. The rationale for 
making such a hypothesis is presented throughout the remain- 
der of the paper. 
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Fig. 2. Montage of VUV spectra of B-type supergiants recorded 
with the International Ultraviolet Explorer (Tt/Ej^Spectra 
downloaded from the MAST Scrapbook interactive web site.) 
The spectral region shown here includes the C II /Lil335, 1336 
resonance doublet. Stars in adjacent spectra are relatively close 
in temperature. 

3. Identifying the proposed P Cygni mechanism 
with stimulated resonance Rayleigh scattering 

In any type of stimulated scattering mechanism, one would ex- 
pect both momentum and energy to be conserved in the unit 
step of the process. As already outlined in the Introduction, in 
the unit step of the process here being invoked to explain P 
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Cygni profiles, three events simultaneously occur at the posi- 
tion of an ion being accelerated in the solar wind. Q A con- 
tinuum photon propagating radially outwards from the star, and 
having an arbitrary frequency v t that is blueshifted with respect 
to the rest frame frequency v a of a resonance transition of the 
ion, is absorbed by the latter. A photon of frequency V2 is 
emitted by the ion in the direction pointing radially inwards to- 
wards the star. The radial velocity of the ion is increased 
from v to v + Av. 

Let relativistic corrections initially be neglected. Requiring 
that energy be conserved in the unit scattering step implies that 

/ivi = hv2 + m/v(Av), (1) 

while requiring that momentum be conserved dictates that 

— + — = mj(Av). (2) 
c c 

From these equations, it then follows that 



vi = v 2 



(i-i)" 



(3) 



Equation Q implies that energy and momentum could in prin- 
ciple both be conserved, even if the value of v 2 were to differ 
in each scattering event. However, if the photon emitted at V2 
were always at the same fixed frequency, this would allow the 
V2 light wave to become amplified much more as it propagates 
radially inwards towards the star, since the rate of stimulated 
scattering at any point in the solar wind is proportional to the 
V2 light wave intensity at that same point. The latter would ob- 
viously be much greater if the backwards-emitted photon in 
each unit step were always at the same fixed frequency. In the 
proposed model, not only is this assumed to be the case, but it 
is also postulated that V2 = v„. 

One can also attempt to view the unit step in the proposed 
scattering process simply as spontaneous resonance Rayleigh 
scattering (i.e. elastic scattering) occurring in the rest frame of 
the ion being accelerated. Viewed in this rest frame, the pump- 
ing frequency (i.e. the incident continuum photon frequency) is 
vi - vi(v/c), since the ion is receding from the star at velocity 
v. (Here again, relativistic corrections are being ignored.) Thus, 
in the ion rest frame, an oscillating dipole moment would be 
induced at the same frequency, vi - vi(v/c). This oscillating 
dipole moment would radiate in all directions. Photons that are 
emitted in the backwards direction (i.e. towards the star) would 
therefore be at V2 = [vi - Vi(v/c)] - [vi - Vi(v/c)] (v/c), imply- 
ing that 



vi = 2 • (4) 

[1 - (v/c)] 2 

The values of Vi given by Eqs. Q and (0} are the same to within 
a factor - (v 2 /c 2 )J , which would result in a difference be- 
tween these quantities of only about one wavenumber at typi- 
cal resonance frequencies (~ 100,000 cirT 1 ) of solar wind ions 
moving at ~ 1000 km/sec. Although the effects of this differ- 
ence on the proposed P Cygni scattering model would be very 
small, it is nonetheless of some interest to compare the two sce- 
narios when relativistic effects are taken into account. We here 



consider first the changes that would occur in the spontaneous 
Rayleigh scattering picture. Use of the known expression for 
the relativistic Doppler shift would imply that, in the rest frame 
of an ion in the solar wind moving away from the star with ve- 
locity v, the frequency of an incident continuum photon at V\ 
would be 



V = Vl 



4 



= Vl 



1 + 



(5) 



Therefore, the frequency of photons emitted in the direction of 
the star by a dipole moment oscillating in the ion rest frame at 
the frequency V would be 



1 - v - 



V" 



Interestingly enough, the above equation expresses exactly the 
same relationship between v\ and V2 that is represented in Eq. 

To complete the comparison of scenarios, one should fi- 
nally consider the equations for conservation of energy and 
momentum that would hold in the unit step of the proposed 
scattering process when relativistic effects are included. The 
equation equivalent to Eq. Q is: 



/ v 2 r 3/2 
hv\ = hvi + m/v II -\ (Av), 

while that equivalent to Eq. (0 is: 

2\-3/2 
A 



hv\ hv2 L v" 

1 mAl — 

c c \ 



(Av). 



(6) 



(7) 



From Eqs. l|6) and Q, an equation identical to Eq. (0 follows. 
Thus, when relativistic effects are included, the equations relat- 
ing vi and V2 in the unit step of the proposed scattering process 
and in spontaneous Rayleigh scattering are seen to be identical, 
strongly implying that the former should be identified with the 
latter. As has already been emphasized, in the proposed model 
it is assumed that V2 = v , and that the presence of an intense 
light wave at this same frequency at the position of every ion in 
the solar wind preferentially stimulates (i.e. induces) the pho- 
tons produced by the Rayleigh scattering process to be emitted 
in the backwards direction, while simultaneously causing the 
ion to be accelerated radially outwards. While it is hard to pro- 
vide the conditions necessary to unambiguously demonstrate 
stimulated Rayleigh scattering in the laboratory, it is nonethe- 
less quite apparent that the astrophysical environments of hot 
stars which display P Cygni profiles would offer optimum con- 
ditions for this particular type of stimulated scattering to oc- 
cur. 

In the model, all photons produced by the stimulated 
Rayleigh scattering process throughout the entire volume oc- 
cupied by the solar wind occur at the same monochromatic 
frequency v a and propagate radially inwards towards the star. 
All photons nonlinearly absorbed by the stimulated scattering 
process are continuum photons that are emitted from the pho- 
tosphere of the illuminating star and propagate radially away 
from it. The frequencies of these absorbed photons, which at 
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all points in the solar wind supply the pumping energy needed 
to drive the stimulated Rayleigh scattering process, span a wide 
spectral range that originates at v and extends to higher ener- 
gies. From Eq. an ion moving with radial velocity v will 
nonlinearly absorb continuum light in a narrow frequency band 

centered at v\ = v [\ + (l - - v + 2v (v/c). As an ion 
becomes accelerated more and more, the continuum light it ef- 
fectively absorbs occurs at higher and higher frequencies. One 
thus can see how the stimulated Rayleigh scattering process is 
able to efficiently convert most of the incoherent light contin- 
uum photons emitted from a star over a wide, blueshifted spec- 
tral range into the same number of essentially monochromatic 
coherent light photons in the form of a spherical light wave at 
v a that radially converges upon the star's photosphere. 

From the above paragraph, one can now comprehend the 
basis for the statement made in the Introduction that use of the 
standard relationship v B - v a — v v max /c to determine the max- 
imum ion velocity v max attained in a star's solar wind results in 
an overestimation of this quantity by almost exactly a factor 2. 
However, this statement does assume that the most blueshifted 
absorption occurring in a P Cygni band almost entirely repre- 
sents nonlinear absorption. 

From Eq. |7} one has that 

2kVo (8) 



Av 



cm i 



This equation states that the velocity increase Av occurring in 
each unit scattering event is always roughly the same - that is, 
to a first approximation, Av is not a function of the velocity v of 
the ion involved in the scattering event. It depends only on in- 
herent properties of the ion being accelerated in the solar wind. 
In the case of the C II ion P Cygni profiles shown in the two 
lower spectra of Fig. [2 one has that Av a 50 cm/sec. Although 
each nonlinear scattering event results in only a modest veloc- 
ity increase for the ion being accelerated, the rate of occur- 
rence of such events will be very large in the region of the so- 
lar wind where the particles are being significantly accelerated, 
due to the Rayleigh scattering process becoming stimulated. In 
the case of Fig. [2;, the short wavelength limit of the P Cygni 
absorption region appears offset from the shortest-wavelength 
C II doublet component by about 200 cm -1 . Thus, according 
to the proposed nonlinear model, v max would here be about 395 
km/sec. An ion accelerated to the maximum velocity v max in the 
solar wind would therefore have had to participate in at least ~ 
790,000 unit step scatterings. 

4. Outline of a model for stimulated Rayleigh 
scattering occurring in P Cygni stars 

In principle, one should be able to model the proposed stim- 
ulated scattering scenario with a set of differential equations 
involving a number of dependent and independent variables. 
Equations for a model possessing minimum mathematical 
complexity are outlined below, and some discussion is given 
of the parameters involved. However, at this stage only qualita- 
tive statements regarding this model can be made, as attempts 
have not yet been made to provide computer-based numerical 
solutions for even this rudimentary system of equations. 



In the conceptually simplest type of model, there would be 
only one independent variable. This would be r, the radial dis- 
tance from the center of the illuminating star. One primary de- 
pendent variable would be v(r), the velocity of an ion in the 
solar wind at radial position r. Assuming that there is only a 
single ion velocity associated with each value of r effectively 
presupposes that the nonlinear ion acceleration process {i.e. 
stimulated Rayleigh scattering) commences at some given ra- 
dius R > R (R being the star's photospheric radius), and that 
at this radius a spherically uniform, radially expanding flow 
of ions occurs, with each ion in the flow moving at the same 
radial velocity v(R ). Let K2 represent the total rate at which 
ions of a given P Cygni species continually cross the surface of 
the sphere at R„, i.e. K2 = AnR^ni{R )v{R ). Here nj(R ) is the 
value at R of another primary dependent variable, the ion den- 
sity at r, «/(r). Steady-state conditions are assumed to prevail in 
the model, and since it is also postulated that all of the ions es- 
cape into interstellar space, an effective equation of continuity 
must exist. One therefore would have 



47rr 2 n/(r)v(r) = K 2 . 



(9) 



The third primary dependent variable is faif), the flux at r of 
the monochromatic, radially inwardly propagating laser wave 
at V2 = v a . The growth of fair) occurs due to stimulated 
Rayleigh scattering, and can be represented by the equation 



dr 



oir)ni{r)<p\{r)4>z{r). 



(10) 



With this equation, two additional dependent variables, <pi(r) 
and cr(r), appear to have been introduced. However, it is 
straightforward to write down the functional dependence on 
r of one of these variables, and the other variable can easily 
be expressed in terms of the previously introduced dependent 
variable v(r). Both these variables will now be discussed. 

The function (f>\(r) is the star's continuum flux at position r 
that is effective in pumping the stimulated Rayleigh scattering 
process occurring at that point. Basically, it represents the con- 
tinuum light at r contained within a certain fixed narrow band 
width Av„. This band width Av„ should approximately corre- 
spond to the spectral line width of the inwardly propagating 
V2 laser radiation. Initially, it can be regarded as an adjustable 
parameter. On the basis of the proposed model, one can write 
(pi (r) = Ki/r 2 , with K\ being a constant easily determinable 
from the emissive properties of the star generating the solar 
wind. The reason why it is possible to write <p\(r) as a sim- 
ple inverse square function of r follows from the discussion 
given in Sec. III. At any radius r, the only pump light that is 
effectively absorbed occurs in a narrow band width Av„ that is 
spectrally located at a v a + 2v (v/c). Thus, nowhere in the solar 
wind is the continuum pump light depleted by ions located at 
positions with smaller r values. 

The function cr(r) is an effective two-photon cross-section 
for the stimulated Rayleigh process. Its functional dependence 
on r is given by 



o-(r) 



(vi(r)[l-*2]-v f 
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u„c^ 



(11) 



y2y 2 (r)' 

this equation simply reflecting the fact that excitation is occur- 
ring on the wing of a lifetime broadened Lorentzian resonance 
profile. The constant cr e can be determined from the radiative 
properties of the P Cygni resonance transition. 

A third independent equation relates two of the three inher- 
ently dependent variables to each other. The radial velocity v(r) 
of each ion located in the volume between r and r + dr will be 
increased by an amount cr(r)(pi(r)(p2(r)(Av) per second. Since 
it takes a time dr/v(r) for an ion to reach r + dr, one can write 



v(r) 



dv(r) 
dr 



cr(r)Mr)Mr)(M. 



(12) 



From Eqs. ij^Ji, 1 11 Oi l, and dl2l >. it appears that obtaining a so- 
lution to the system of equations here discussed amounts to 
solving the following two coupled equations: 



2(r) 



dr 



<h(r) 
rV(r) 



and 



vV) 



dv(r) 
dr 



{hv)(r (?Ki 



<h{r) 



(13) 



(14) 



All the constants appearing in brackets in these two equations 
are in principle known for transitions having P Cygni line 
shapes in stars such as £ Pup. It would be interesting to see 
if numerical solutions can indeed reveal a threshold type of be- 
havior, with parameter values in certain ranges resulting in the 
occurrence of "lasing" in the star's solar wind. 

5. Possible spectroscopic consequences of the 
proposed model 

5.1. The fluorescence component of a P Cygni profile 

As has been indicated above, a signature feature of the pro- 
posed nonlinear P Cygni model is an intense monochromatic V2 
laser wave that radially impinges upon the star's photosphere. 
Since this wave would be directly unobservable in any line of 
sight, one should therefore consider what kinds of observable 
secondary effects might result from its presence. 

One predictable effect would be fluorescence seen in the 
line-of-sight that is emitted by P Cygni ions moving in the 
solar wind. As discussed in various texts (e.g. Loudon 2000 1, 
when there is no collision broadening, and the excitation inten- 
sity is weak, the Lorentzian cross-section for linear resonant 
light scattering by a two-level atom exclusively applies to the 
elastic (i.e. Rayleigh) component. The probability of inelastic 
scattering is zero, i.e. no real linear excitation of the atom oc- 
curs. The situation becomes entirely different, however, when 
intense excitation (as produced by a laser beam, for example) 
is present. Then real excitation of the upper level of the atom 
can result, even when the exciting light is non resonant (c.f. 



Knight & Milonni 1980 1. One thus could expect ions in the so- 
lar wind to undergo real transitions to the upper level of the 
P Cygni resonance, driven by an intense inwardly propagat- 
ing laser wave at v%. Such ions would then fluoresce. At first 
glance, it might appear that the fluorescence occurring at v > v G 
would be absorbed in pumping the stimulated Rayleigh scatter- 
ing processes that occur in the solar wind. However, it can be 
readily shown that this cannot be the case. Consider an ion in 
the solar wind moving with velocity v in a direction towards 
Earth. If this ion becomes inelastically excited by the V2 laser 
wave, the fluorescence it immediately would radiate towards 
Earth would occur at the frequency v„ + v (v/c). However, it 
was shown in Sec. Ill that the continuum pump light necessary 
to make the same ion undergo stimulated Rayleigh scattering 
occurs at frequency v + 2v (v/c). Thus, the fluorescence emit- 
ted by the ion towards Earth cannot be nonlinearly absorbed 
anywhere in the solar wind. (In principle, it could be linearly 
absorbed to some extent by surrounding ions in the solar wind, 
but in the proposed nonlinear model this process is being en- 
tirely neglected.) Fluorescence in the line-of-sight occurring at 
v <v will neither be linearly nor nonlinearly absorbed. 

In view of the deduction made in the last paragraph, it 
would appear that, on the basis of the proposed nonlinear 
model, a rough replica of a typical P Cygni profile could 
be conceptually constructed in the following manner. Start 
with the basic continuum level of the star. Subtract from this 
the blueshifted light that is absorbed in pumping stimulated 
Rayleigh scattering. To this depleted continuum then simply 
add the fluorescence intensity radiated towards Earth by all 
the ions in the solar wind. Such fluorescence will generally be 
much stronger for v > v a than for v < v a , due to occultation by 
the star. Via this conceptual construction, a simple explanation 
is therefore provided for a characteristic observed feature of 
P Cygni profiles, namely, that the strongest attenuation of the 
continuum appears to occur at frequencies that are significantly 
offset from v„. In reality, strong attenuation of the continuum 
may also be occurring much closer to v , but this is masked by 
the presence of a strong fluorescence signal. 

From the fluorescence mechanism here proposed, it also 
follows that the reddest fluorescent emission observed should 
have an absolute value frequency offset from v no greater than 
half that of the bluest P Cygni absorption. In most P Cygni 
spectra this seems approximately to be the case. However, it is 
perhaps well to note here that fluorescence appears not to be as 
intrinsic a feature of P Cygni profiles as is extended blueshifted 
absorption. For the P Cygni spectra observed in Morton ( 1976 1, 
it was noted that the equivalent width of the emission com- 
ponent is always significantly less than that of the absorption 
component. In the P Cygni profile of N III in ( Pup, in fact, 
the fluorescence component is completely absent. It was sug- 
gested in Morton dl976> that this might be the result of absorp- 
tion occurring close enough to the stellar surface to hide the 
emission via occultation by the star. To obtain even a rough 
estimate of the expected fluorescence profile based upon the 
nonlinear model would require at least semi-quantitative solu- 
tions of equations such as those discussed in Sect. 4 to be made. 
From the present section, it would appear that one should mod- 
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ify those equations by including a loss term for the variable 
<j>2(r), in order to represent strong fluorescence being excited. 

5.2. Two-photon absorption bands 

In principle, one could seek to directly substantiate the pro- 
posed nonlinear model by clearly identifying in a P Cygni spec- 
trum Doppler broadened absorption bands representing res- 
onantly enhanced two-photon absorption of continuum light 
from the star by various ion species present in its photosphere, 
or just outside. Such two-photon absorptions could in princi- 
ple be induced if powerful, monochromatic beams at P Cygni 
ion rest frame frequencies v were indeed incident upon the 
star's photosphere, as the nonlinear model presupposes. For 
such two-photon absorptions to be strong enough to be de- 
tectable, near resonances would have to exist between some of 
the P Cygni rest frame frequencies and some of the transitions 
of the two-photon-absorbing ions. 

The authors have attempted to identify such two-photon ab- 
sorption bands in the case of the P Cygni star £ Pup. With use 
of the NIST Atomic Spectra Database web site, spectral inter- 
vals extending ~ 200 cirT 1 on either side of the rest frame fre- 
quencies v of the strongest P Cygni profiles shown in Morton 
J1976I were examined for likely absorbing-ion candidates, but 
no convincing ones were found. The spectral intervals exam- 
ined were as follows: C III (977 ± 2A); C IV (1551 + 5 A, 1548 
+ 5A); N V (1243 ± 3.2A, 1239 + 3.2A); Si IV (1403 + 4A, 
1394 ± 4A); N III (992 + 3 A, 990 + 3A); S VI (945 + 5 A, 933 
+ 5A); and O VI (1032 + 5A, 1038 ± 5A). The most promis- 
ing coincidence found was between the wavelength (1038 A) 
of one of the O VI doublet components and a strongly allowed 
C II ion transition originating from a level only 63 cirT 1 above 
its ground state (Fig. [3}- The difference between the O VI reso- 
nance and C II transition frequency is in this case only 55 cm -1 . 
With this scheme (and with another somewhat less resonant 
one also shown in Fig. O, one predicts two-photon absorption 
bands to exist at 1381.7 A, 1382.08 A, 1382.94 A, and 1383.3 
A, with the latter two being stronger than the former two. 

On the basis of several high resolution VUV Copernicus 
scans of ( Pup recorded in 1973-1975, interstellar absorption 
lines seen in the line-of-sight to this star were identified and 
tabulated in Morton (1978 1. This same publication also con- 
tains a table of 52 lines towards ( Pup that were unidentifi- 
able in 1978. In 1991, Morton published a large compendium 
(Morton fl 99 It of interstellar absorption lines based upon data 
obtained from many space objects with use of the Hubble 
Space Telescope. Included in this finding list was an updated 
version of the original unidentified lines list of Morton ( 1978 1, 
the number of such lines towards f Pup by 1991 having been 
reduced to 42, largely as a result of the discovery of seven new 
Fe II lines that were not known in 1978. Still remaining unas- 
signed in 1991 were four lines in £ Pup occurring as an iso- 
lated group at 1382.305 A (8.9), 1382.593 A (22.1), 1383.205 
A (30.8), and 1383.720 A (18.6) (equivalent widths in mA in 
parentheses). One sees that the four unidentified lines in ( Pup 
are tantalizingly close to the four predicted two-photon absorp- 
tion bands. However, despite the existence of these intriguingly 
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Fig. 3. C II partial energy level diagram showing two-photon 
absorption transitions that are likely to be induced by appli- 
cation of strong monochromatic coherent radiation at 96375 
cm -1 , the frequency of the O VI resonance at 1038 A. 



close wavelength matches, it would be apparently difficult to 
argue that C II would not be fully ionized in the photosphere of 
a star as hot as ( Pup. Strong C II absorptions are seen in the 
£ Pup spectrum, but astronomers believe that these represent 
interstellar lines. 

It is of interest here to comprehend roughly what the max- 
imum intensity of the monochromatic V2 wave could be as 
it impinges upon a star's photosphere. We here consider the 
case where V2 corresponds to the O VI P Cygni resonance at 
1032 A in ( Pup. In this star, the corresponding blueshifted 
absorption region extends to about 1020 A. Under the assump- 
tion that roughly 80% of the continuum photons emitted by 
the star between 1020 A and 1032 A are nonlinearly absorbed 
and entirely converted to photons of the inwardly propagating 
monochromatic wave at V2, it follows that the intensity of the 
latter would be roughly 200 kW/cm 2 as it impinges upon the 
star's photosphere. In this estimate it is assumed that the tem- 
perature of £ Pup is 50,000K, and that this star emits as a per- 
fect blackbody. It is also here assumed that no loss of V2 wave 
intensity occurs as a result of fluorescence excitation. As out- 
lined in Sect. 5a, one actually expects such loss to be quite 
large. 

6. Linear vs. nonlinear P Cygni mechanisms in 
luminous stars 

As noted above, current theoretical models for line driven 
winds of luminous hot stars are all based upon the effects of 
linear resonant absorption of continuum light emitted by such 
stars. In such models, the essential physical processes that oc- 
cur are in many ways similar to those that characterize the 
nonlinear model described in Sect. 3. (Chapter 8 of Lamers & 



8 



P. P. Sorokin & J. H. Glownia: A model for selective line-driven acceleration of ions 



Cassinelli ( 1999 1 is a good exposition of the currently accepted 
linear theory of line driven winds.) For example, it is easy to 
demonstrate that the outwardly directed momentum transfer to 
an atom or ion in the solar wind in a linear scattering event is 
liVg/c - i.e. half the value occurring in the unit step of the non- 
linear process. In the linear theory, just as in the nonlinear case, 
the Doppler shifts of atoms or ions moving in the outer portions 
of a solar wind in principle allow the atoms to absorb undimin- 
ished continuum photons in their line transitions. However, the 
light that effectively drives the acceleration of an atom mov- 
ing with a given radial velocity v in the case of the nonlinear 
mechanism is blueshifted from v twice as much as in the lin- 
ear case. This implies that occurrence of the former mecha- 
nism {i.e. stimulated Rayleigh scattering) would have already 
strongly depleted the supply of continuum photons necessary 
to drive the linear process everywhere in the solar wind. Since 
the stimulated Rayleigh scattering process should occur at its 
highest rate very close to the photosphere of the illuminating 
star, it is even possible that in very bright stars it completely 
prevents linear scattering from being an effective acceleration 
mechanism. 
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